Draft version September 29, 201 1 

Preprint typeset using I4TgX style emulateapj v. 08/22/09 



NONLINEAR SMALL-SCALE DYNAMOS AT LOW MAGNETIC PRANDTL NUMBERS 



o 

(N 

D 

00 
(N 

I 

o 



(N 
> 

o 



Axel Brandenburg 

NORDITA, AlbaNova University Center, Roslagstullsbacken 23, SE-10691 Stockliolm, Sweden; and 
Department of Astronomy, Stockholm University, SE-10691 Stockholm, Sweden 
Draft version September 29, 2011 

ABSTRACT 

Saturated small-scale dynamo solutions driven by isotropic non-helical turbulence are presented at low mag- 
netic Prandtl numbers Pim down to O.OL For Ptm < 0.1, most of the energy is dissipated via Joule heat and, in 
agreement with earlier results for helical large-scale dynamos, kinetic energy dissipation is shown to diminish 

1 /2 

proportional to Pr^^^ down to values of 0. L In agreement with earlier work, there is, in addition to a short 
Golitsyn fc^^^/"^ spectrum near the resistive scale also some evidence for a short k^^ spectrum on larger scales. 
The rms magnetic field strength of the small-scale dynamo is found to depend only weakly on the value of Ptm 
and decreases by about a factor of 2 as Pim is decreased from 1 to O.OL The possibility of dynamo action at 
PrM = 0. 1 in the nonlinear regime is argued to be a consequence of a suppression of the bottleneck seen in the 
kinetic energy spectrum in the absence of a dynamo and, more generally, a suppression of kinetic energy near 
the dissipation wavenumber. 
Subject headings: MHD - turbulence 



1. INTRODUCTION 

In astrophysical turbulence, dissipation of kinetic and mag- 
netic energies tends to occur on length scales much shorter 
than the scale of the energy-carrying eddies. Even though 
both kinetic and magnetic dissipation scales are compara- 
tively short, the current indications are that it does matter 
which of the two is the shorter one and by how much. Their 
ratio is the magnetic Prandtl number, Ptm . For stars and liq- 
uid metals we have Pim ^ L while for galaxies Pim ^ 1. 
An important example where the value of Pym is believed to 
matter is the small-scale dynamo that converts kinetic turbu- 
lent energy into magnetic energy under isotropic conditions. 

A dynamo is only possible when the energy conversion is 
efficient and larger than the magnetic energy dissipation. This 
is quantified by the magnetic Reynolds number, Rcm, which 
is a nondimensional measure of the inverse magnetic dissipa- 
tion rate. The critical value of Rsm, above which dynamo 
action occurs, is known to increase with decreasing values of 
PrM feogachevskii & Kleeorin 1997; Boldyrev & Cattaneo 
2004; Haugen et al. 2004f ISchekochihin et alj 12004 120051 
2007; Iskakov et al. 2007^. In the following we define the 
magnetic Reynolds number as Rejv/ = Urms/??^/, where 
Urms is the rms velocity fluctuation of the turbulence, rj is 
the magnetic diffusivity, and kf is the wavenumber of the 
energy-carrying eddies, i.e., the wavenumber where energy 
is injected into the system. The critical value of Re a/ is 
then found to be around 35 for P rM = 1 and arou n d 100 
for PrM = 0.2, but note that for iSchekochihin etaTI (12004 
I2005L 120071) and llskakovetalJ (|2007|) the values of ReM are 
defined such that they are about 1.5 times larger than those 
used here or in Haugen et al. (2004). For PrM = 0.1, how- 
ever, no small-scale dynamo action has yet been found. This 
may easily be a limitation of not having been able to increase 
the fluid Reynolds number. Re = Re m /PrM, beyond 2000 , 
which limi ts ReM to 200 for Ptm = 0.1 (llskakov et al.ll2007l: 
ISchekochi hin et al■l2007^ . Larger values of Re have been pos- 
sible by using hyperviscosity, giving access to larger values of 
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Re and smaller values of P im for fixed Re^- In that case 



llskakov et all (120071) and ISchekochihin et alJ (l2007h found 
small-scale dynamo action for PrM = 0.05 and ReM = 150, 
i.e., the dynamo is now easier to excite than for PrM = 0.1. 
The reason for this is believed to be connected with the fact 
that the properties of small-scale dynamos depend on the ki- 
netic energy spectrum at the resistive scale. For PrM = 1, this 
scale is the viscous scale where the velocity field is smooth in 
the sense that the velocity difference 5u over a separation 5£ 
scales linear ly, i.e., 5u ~ 6£. For Pr]\i ^ 1, following the 
argument of iBoldvrev & Cattaneo! (2004), the resistive scale 
falls in the inertial range where the velocity field is rough and 
Su ^ 61'' with C « 0.4, so the velocity field would not be 
differentiable, making dynamo action inefficient. However, 
for PiM = 0.1, the kinetic energy spectrum is even shal- 
lower than in the inertial range, so the local value of C is 
even smaller and the velocity field rougher than in the iner- 
tial r ange. This phenomenon is kriown as the bot tleneck ef- 
fect (lFalkovichll994l;lKaneda et~ani2003l;lDobiere t al. 2003). 
This bottleneck effect is believed to be the reason why i?m,crit 
reaches a maximum at PrM ~ 0.1. 

The usage of hyperviscosity does exaggerate the bottleneck, 
which still exists even for the regular viscosity operator. It 
would therefore be useful to verify small-scale dynamo ac- 
tion for small values of Prj\/ using the regular viscosity op- 
erator. This will be done in the present paper. In addition, 
we shall consider here the nonlinear regime, which has the 
advantage that at small values of PrM, much of the kinetic en- 
ergy is diverted to magnetic energy before it is dissipated vis- 
cously. This allows one to increase the fluid Reynolds num- 
ber beyond the maximal value that would normally be possi- 
ble at a given resolution. This has been demonstrated in the 
conte xt of helicity-driven large-scale dynamos (Brandenburg 
i2009h . whose onset conditions are essentially independent of 
flie value of PrM (Br andenburg 200 L 2009; Mininni2007,) . 
This is not the case for the small-scale dynamos considered 
here, where the flow is statistically isotropic and non-helical. 
Our strategy for reaching low values of Pr^f is the same as 



ng 

that of iBrandenburgj (l2009t) . We start with a simulation of a 
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Fig. 1 . — Visualizations of and Uz for PrM =0.01 (left), Pfm = 0.02, and Ptm = 0.05. All rans are for Re a/ 160 using 512^ mesh points. 



saturated small-scale dynamo at Ptm = 1 and then increase 
the value of Re while keeping the value of Rsm in the range 
150-160, provided the dynamo is still excited. We are here 
particularly interested in the dependence of the saturation field 
strength on Pim and the dissipation rate. 

2. THE MODEL 

Our model is similar to that presented in iBrandenburd 
( l200lll2009h and lHaugen etaTI (I200l l200?), where we solve 
the hydromagnetic equations for velocity U, density p, and 
magnetic vector potential A, in the presence of an externally 
imposed non-helical forcing function /, for an isothermal gas 
with constant sound speed Cg, i.e.. 



dU 

— = -U-VU-clV\np+f+{JxB+V-2pvS)/p, (1) 



9 In p 
dt 



= -U-Vlnp-V -U, 



dA 
'dt 



= U X B rjfioJ. 



(2) 



(3) 



Here, S^- = \{Ui_j + Uj^t) — ^Sij'V ■ U is the traceless rate 
of strain tensor, i/ is the kinematic viscosity, B = "V x A 
is the magnetic field, J = V x B/po is the current den- 
sity, and po is the vacuum permeability. We consider a triply- 
periodic domain of size L^, so the smallest wavenumber in 
the domain is ki = 2tt/L. The forcing function consists of 
plane waves with wavevectors fc whose lengths lie in the range 
1 < |fe|/fci < 2 with an average of fc/ w 1.5 fci. The ampli- 
tude of / is such that the Mach number is Unns/cg ~ 0.1, so 
compressive effects are negUgible (.Dobler et al...2003i) . 



Unless a simulation has been restarted from a previous one 
at another value of Pr^/, we start with a weak Gaussian dis- 
tributed field in all three components of A, zero initial ve- 
locity, and uniform initial density, p — po = const, so the 
volume-averaged density remains constant, i.e., (p) = po- 

In our simulations we vary the fluid Reynolds number and 
the magnetic Prandtl number. 



Re = Urms/j^fc/, PrM = t^/?7, 



(4) 



such that Rejv/ = u^ms/'H^f is in the range 150-160. We 
also present a few results for Rcm around 220. We monitor 
the resulting kinetic and magnetic energy dissipation rates per 
unit volume. 



eK = (2i/pS^), EM = (r;/ioJ^), 



(5) 



whose sum, er — + ^m, is the total dissipation rate. We 
use the fully compressible PENCIL Code' for all our calcu- 
lations. We recall that, for the periodic boundary conditions 
under consideration, {23'^) = {W^) + |((V • t/)^), high- 
lighting thus the analogy between vorticity W — 'V x U and 
J in the incompressible and weakly compressible cases. 

3. RESULTS 

3.1. Small-scale magnetic and velocity features at low Pim 

In Figure [U we present visualizations of and Uz on the 
periphery of the domain for three runs with Pim ranging form 
0.01 to 0.05. Even though the value of ReM is the same in 
all three runs, the magnetic field seems to have smaller scale 
structures in the low Pr^/ case. The appearance of smaller 

' |http://www.pencil-code.googlecode.co^ 
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scale structures is particularly clear in the visualization of the 
velocity field for Ptm = 0.01. 

In Table [T] we summarize some essential properties of 
the simulations for a sequence of simulations with differ- 
ent values of Pr^ between 0.01 and 1, but similar val- 
ues of Rcm of around 150-160. The rms field strength 
relative to the equipartition value, i?cq — Urms-^/PoMo. is 
about 0.3 for 0.05 < Pvm < 0.2, while for the runs with 
PrM = 0.02 and 0.01 it is about 0.15 and 0.12, respec- 
tively. This is still a remarkably weak dependence that was 
not ex pecte d based on the earl ier res ults by Iskakov et al] 
( 120071) and ISchekochihin et all ( 120071) for the onset condi- 
tions of the small-scale dynamo. As Ptm is decreased from 
1 to 0.01, ex decreases and cm increases. However, the runs 
for Rej\/ ~ 160 are rather close to the onset of dynamo ac- 
tion. This becomes clear when comparing with two other 
runs for Rej\/ = 220 and Pr]\i = 0.1 and 0.02; see Ta- 
ble m For Pvm = 0.1, Sims/^oq ~ 0.32 and the ratio 
eK/er has dropped from 0.39 to 0.24, while for Pvm = 0.02, 
Binis/Bcq ~ 0.34 and the ratio eK /^t has dropped from 0.6 
to 0.08. Thus, we see that for values of Re7\/ that are not too 
close to the onset of dynamo action, the Pr^/ dependence of 
Brms/Bcq is negligible and ex continues to drop. 

The magnetic dissipation wavenumber, /cm = {^m /v^Y^^^ 
is about 30 for all runs with Rcm ~ 150-160, while the ki- 
netic dissipation wavenumber, fc/f = {ck /i^^Y^^, increases 
gradually with decreasing values of Prj\/ (or increasing val- 
ues of Re). 

3.2. Spectral properties and energy dissipation 

We consider here kinetic and magnetic energy spectra, 
Exik) and EM{k), respectively. They are normalized in the 
usual way such that / Ek dk ~ ^po{U'^) and J EM<ik = 
^Pq^{B'^). In Figure |2] these spectra are compensated with 

ey^/^fc^/^ For Ptm = 0.02 and 0.01, the kinetic en- 
ergy spectra show a clear bottleneck effect, i.e., there is a 
weak upris e of the compensated spectra toward the dissipative 
subra nge (IFalkovichlll994l: iKanedaet al]l2003t iDobler etalJ 
l2003h . The compensated magnetic energy spectra peak 
around k — 20/ci . Both toward larger and smaller values of 
k there is no clear power law behavior, although the slopes of 
the spectrum of iGohtsvnl ([i960): Moffatt (196 j]) and 

the scale-invariant k~^ spectrum ( | Ru zmaikin & Shukurov' 
Il982t iKleeorin & Rogachevskiilll994n Kleeorin et al. 1996) 
are shown for comparison. 

It turns out that for small values of Pim, dynamo action is 
maintained for Rcm ~ 160, corresponding to Re w 7800. 
This value of Re is rather large for a resolution of 512'^ mesh 
points and one must be concerned about insufficient resolu- 
tion. Similar circumstances were encountered previously in 
connection with simulations of large-scale dynamos at low 
values o f PrM ([B randenburg 2009), and even at large values 
of Pim (iBrandenbur g 201 1). In the former case, much of the 
energy dissipation occurs magnetically via Joule dissipation, 
leaving thus very little energy in the rest of the kinetic en- 
ergy cascade. This allows us then to decrease i/ further, while 
still allowing the remaining kinetic energy to get dissipated. 
However, kinetic and magnetic energies are quite intermittent 
(uppermost panel of Figure |3]l and there can be extended pe- 
riods over which the magnetic energy drops well below the 
kinetic energy. Nevertheless, the magnetic energy dissipation 
is still in excess of the kinetic energy dissipation; see Figure[3] 

For PrM — 0.01, the nominal value of Re is 16,000. 



The kinetic energy spectrum extends now further to higher 
wavenumbers, but it shows still a monotonic decrease down 
to the Nyquist wavenumber at fc = 256fci. As can be seen 
from Table [T] the nominal dissipation wavenumber, kx, is 
now well outside the range of resolved wavenumbers, so it 
is clear that higher resolution would be needed to resolve the 
smallest scales properly. However, as far as the dynamo is 
concerned, most of the magnetic field generation occurs at 
wavenumbers below 30fci, which is where the compensated 
magnetic energy spectrum begins to show a clear decline into 
the dissipation subrange. Until that wavenumber, significant 
amounts of kinetic energy are being channeled into magnetic 
energy via the dynamo, which lowers the kinetic energy dis- 
sipation and is the main reason for being able to run such low 
PrM cases. 

The velocity field is relatively steady over the course of 
the simulation, but the magnetic field and also the magnetic 
energy dissipation vary significantly; see Figure H) How- 
ever, although there can occasionally be a dramatic decline 
in the magnetic field, it tends to recover subsequently, sug- 
gesting that dynamo action is still possible at small values of 
PrM. Obviously, in addition to longer run times, it is nec- 
essary to perform simulations at higher resolution, which is 
not currently feasible if one wants to cover sufficiently many 
turnover times. 

We have already argued that the somewhat erratic behav- 
ior of the dynamo at Re a/ = 160 is a consequence of being 
close to the marginal value. The time evolution for the case 
with Rcm — 220 and PrM = 0.02 is shown in Figure |5] 
Both Bj-yns and ex are now much closer to being statistically 
steady. Furthermore, the value of i?r,ns/^cq is now « 0.3 
both for PrM = 0.1 and for Pi'm = 0.02. This suggests that 
the saturation level of the dynamo is now beginning to be in- 
dependent of the value of Pim ■ 

Earlier work on large-scale dynamos from helical isotropic 
turbulence showed that the total magnetic energy dissipa- 
tion is larger than in hydrodynamic turbulence. This is best 
demonstrated by considering the conventionally defined di- 
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Fig. 2. — Compensated kinetic and magnetic energy spectra for runs with 
PrM = 0.05, PrM = 0.02, and Prjv/ = 0.01 for Re^/ 150 as well as one 
run with Pr^/ = 0.02 and Rbm ~ 220. The resolution is in all cases 512'' 
mesh points. The two short straight lines give, for comparison, the slopes 2/3 
(corresponding to a spectrum for k < 20fci) and —2 (corresponding to 
a A;"^^'''' spectrum for k > 20fci). 
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TABLE 1 

Summary of runs for different values of Prjj/ and Re^/ 



150-160. 



PrAf 


ReM 


-Brms / -Bcq 


a 






kK 




At/r 


Res. 


0.01 
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0.12 ±0.02 


0.34 ±0.03 


0.49 ±0.13 


0.68 ±0.06 
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Fig. 3. — Root-mean-square velocity, ratio of magnetic to kinetic energy, 
as well as the ratio of magnetic to kinetic energy dissipation for the run with 
Prjv/ = 0.02 using 512^ mesh points. In the last two panels the dashed lines 
denote normalization with respect to the instantaneous values of Boq and e^, 
respectively, while the solid lines refer to normalizations based on the time 
averaged values of B^q and ejf. 

mensionless dissipation parameter 

where U is the one-dimensional rms velocity, which is related 
to Wrms via — u^j,jg/3, and L is the integral scale which is 
related to kf via jir/kf. In non-helical turbulence this value 
is typically around 0.5 (see, e.g. Pearson et al. 2004), but in 
helical turbulence with large-scale dynamo action this value 
is around 1.4; see also Brandenburg (2011). In Figure |6] we 
use time averaged dissipation rates, which, for simplicity, are 
also denoted by £k, em, and er- In the upper panel of Fig- 
ure|6]we show that, in the present case of small-scale dynamo 
action from non-helical isotropic turbulence, this value is now 

TABLE 2 

Comparison of runs for Rejv/ ~ 220 and two values of Vtm . 



0.02 0.34 ±0.02 0.70 ±0.04 0.08 ± 0.01 0.92 ± 0.01 
0.10 0.32 ±0.03 0.66 ±0.08 0.24 ± 0.03 0.81 ± 0.02 



Fig. 4. — Same as Figure[3] but for the run with Prj\/ = 0.01. 

closer to the hydrodynamic value and is slightly above 0.6 for 
Pi" A/ — 0.02 and Re a/ = 220; see upper panel of Figure |6l 
In the lower panel of Figure|6]we see that the ratio e^f /cm is 
compatible with a Pr°;f dependence, as was found earlier for 
helical hydromagnetic turbulence ( Bran denburg 2009, 201 1). 
However, for Re a/ = 160, ^kI^m levels off at a constant 
value of w 0.4, which is probably an artifact of Rcm being 
too close to the onset of dynamo action. 

3.3. Possibility of subcritical dynamo action 

We recall that we have used here the strategy of generating 
low-Prjvf solutions by gradually decreasing v, and hence in- 
creasing the value of Re. As in the case of helical dynamos 
jBrande nburg 2009), the fact that a turbulent self-consistently 
generated magnetic field is present helps reaching these low- 
Prjv/ solutions. However, the presence of the magnetic field 
also modifies the kinetic energy spectrum and makes it de- 
cline slightly more steeply than in the absence of a magnetic 
field; see Figure|2l This suggests that the velocity field would 
be less rough than in the corresponding case without mag- 
netic fi elds. Following the reasoning of Boldyrev & Cattaneo| 
( 120041) . this should make the dynamo more easily excited than 
in the kinematic case with an infinitesimally weak magnetic 
field. In other words, there is the possibility of a subcritical 
bifurcation where the dynamo requires a significantly larger 
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Fig. 5. — Similar to Figure|3] but for the run witli Rejv/ = 220, using still 
Vvm = 0.02. 
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Fig. 6. — PrM dependence of the dimensionless dissipation rate, Ct, and 
the kinetic to magnetic energy dissipation ratio, ^k/^m- 

value of Pr7\/ to bifurcate from the trivial B = Q solution 
than the value needed to sustain a saturated dynamo. 

In order to check this hypothesis, we perform an experiment 
where the simulation is continued after having down-scaled 
the magnetic field by a factor of 10. The result is shown in 
Figure |7] where we compare the original simulation with the 
one restarted with a 10 times lower field. One sees a grad- 
ual decline of the magnetic field after a brief initial increase 
of the magnetic field. This initial increase is a consequence 
of the reduced feedback from the Lorentz force, allowing the 
velocity to increase slightly above the previous value (see the 
upper dotted line in Figure |7]l. During the next 100 turnover 
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Fig. 7. — Evolution of magnetic and kinetic energies in the main run (solid 
lines) and after rescaling the magnetic field by a factor of one tenth (dotted 
lines). In that case, the resulting decay rate is 0.006urmsfc/. 

times, the decay rate is about O.OOGurmsfc/, which is about 4 
times smaller than the growth rate of 0.025Mrni sfef for a non- 



helica l dynamo at Pim ~ 1 and Re w 150 dHaugen et alj 
120041) . However, the field still seems to recover and shows in 
the end a behavior comparable to that without rescaling. This 
may suggest that at this value of Prj\/ the dynamo may not be 
subcritical after all. 

The possibility of subcritical dynamo action is well known 
in the geodynamo context, where the flo\y is dr iven by ther- 
mal or composition al convection ( Robert sifl 9881) . and for Ke- 
plerian shear flows (iRincon et al.l 20081) . Also in the context 
of dynamos from forced Taylor-Green flows the possibility 
of subcritical dynamos is well known (Ponty et al. 2007). In 
the present context, subcriticality is likely to be linked to the 
steeper kinetic energy spectrum in the low-Prjv/ regime. How- 
ever, because of extended transients, the results for Ptm — 
0.02 shown in Figure|7]remain inconclusive. For Ptm = 0.1, 
on the other hand, Schekochihin et al. (2007) have not seen 
dynamo action in the linear regime when Rej\/ — 160. 

4. CONCLUSIONS 

In the pre s ent p aper we have extended the w ork of 
llskakov et all (l200'7h and ISchekochihin et all (120071) to the 
nonlinear regime of a saturated dynamo. However, while 
in the former (linear) case the dynamo shows signs of a de- 
pression in the range 0.1 < Ptm < 0.2, the nonlinear sat- 
urated dynamo is found to operate nearly unimpededly in 
th e range 0.02 ^ PrM < 1. Furthermore, unl i ke the work 
of llskakovetm (120071) and ISchekochihin et all (l2007h . who 
used hyperviscosity, we have here used regular viscosity with 
the usual Laplacian diffusion operator. As in earlier work on 
helical large-scale dynamos ( Brandenburg .2009). it is possi- 
ble to reach the regime of low PrM by restarting the simula- 
tions from another one at a larger value of Pr7\/, which reduces 
the kinetic dissipation rate proportional to the square root of 
PrM- Furthermore, in contrast to helical large-scale dynamos, 
where dynamo onset is possible for values of Re a/ of the order 
of unity and independently of the value of PrM, we have here 
the situation where the critical value of Re a/ may be larger 
than the value required to sustain the dynamo once it has sat- 
urated. This means that the dynamo could be subcritical and 
might possesses a finite amplitude instabihty at Re m below 
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and around 160. 

1/2 

We note that the Pr^^^ scaling of the kinetic to magnetic 
energy dissipation ratio, /^m, is still not well understood. 
In view of the definitions of €k and eM in Equation (|5]l, it is 
clear that this implies that 



.1/2 



const. 



(7) 



This means that the usually expected scaling for hydrody- 
namic turbulence, v{pW'^) = const, or the hydromagnetic 
scaling r/(/xoJ^) — const, whic h has been con firmed for 
PiM — 1 (see Fig. 8 of Cand elaresi et al.|[201 Ih . is clearly 
not generally valid and needs to be reconsidered. 

Our results for the magnetic energy spectra are consis- 
tent with those of earlier direct numerical simulations by 
ISchekochihin et aLl (120071 ) in that there is a short Golitsyn 
^-11/3 spectrum near the resistive scale, as well as a short 
spectrum on larger scales. Both pro perties have also 
been seen in liquid sodium experiments ( Odier et al.l 119981; 
iBpurgoin et al. 2002) as well as in large eddy simulations 
IPontv et al.ii2004<) . 

The astrophysical relevance of small-scale dynamo action 
is hardly disputed. Even in situations were large-scale dy- 
namo action is possible, like in the Sun, small-scale mag- 
netic fields are seen ubiquitously even in the quiet photo- 
sphere where there is no evidence of any effects from the 
large-scale field (Cattaneo 1999; Vogler & Schiissler 2007; 
iPietarila Graham et al.l 120 lOt iJinet al.1 12011i) . The present 
work now confirms that the small value of the Sun's magnetic 
Prandtl number may not be a problem with this proposal. Al- 
though one may worry that most of the simulations presented 
so far have overestimated the effects of small-scale dynamo 
action by having chosen values of Prj\/ of the order of unity 
(lBrandenburg.,2005.) . it is remarkable that the field strength 



decreases only slightly when we decrease Pi'm from 1 to 0.02, 
provided Rcm is large enough (Re^f ^ 150). Future work 
will hopefully clarify further the relative importance of large- 
scale and small-scale dynamo action in astrophysical bodies 
like the Sun. 

Another aspect that needs to be addressed in future simu- 
lations concerns the magnetic Prandtl number effects on the 
large-scale properties of the turbulence. This concerns in par- 
ticular the turbulent diffusion of large-scale magnetic fields, 
as can be measured by the quasi-kinematic test-field method, 
for example (Brandenburg et al. 2008). Among other things, 
one would like to confirm that the turbulent magnetic diffusiv- 
ity is not affected by the small-scale magnetic field, which is 
a standard result in mea n-field theory dG ruzinov & Diamong 
[T994tlRadIer et al.l2003tlBrandenburg & Subramanianii2005r 
In that case, if Rcm is close to the onset of small-scale dy- 
namo action, one would expect the turbulent diffusion to be 
independent of the value of PrM. However, for large values 
of Re A/, as we have now seen, the effect of Prj\/ on the small- 
scale dynamo is less dramatic. Thus, even if the turbulent 
magnetic diffusivity was affected by the small-scale field, the 
effect could only be weak. 
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